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HIGHLIGHTS 


•  A  chemistry  adaptive  battery  fuel  gauge  is  derived  based  on  the  PDA  algorithm. 

•  The  PDA  algorithm  selects  the  correct  OCV  curve  online  from  a  library. 

•  For  unknown  batteries,  the  PDA-BFG  selects  the  most  similar  OCV  curve. 

•  A  machine  learning  approach  selects  diverse  OCV  curves  to  fill  the  library. 
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This  paper  considers  the  problem  of  state  of  charge  (SOC)  tracking  in  Li-ion  batteries  when  the  battery 
chemistry  is  unknown.  It  is  desirable  for  a  battery  fuel  gauge  (BFG)  to  be  able  to  perform  without  any 
offline  characterization  or  calibration  on  sample  batteries.  All  the  existing  approaches  for  battery  fuel 
gauging  require  at  least  one  set  of  parameters,  a  set  of  open  circuit  voltage  (OCV)  parameters,  that  need 
to  be  estimated  offline.  Further,  a  BFG  with  parameters  from  offline  characterization  will  be  accurate  only 
for  a  “known”  battery  chemistry.  A  more  desirable  BFG  is  one  that  is  accurate  for  “any”  battery  chemistry. 
In  this  paper,  we  show  that  by  storing  finite  sets  of  OCV  parameters  of  possible  batteries,  we  can  derive  a 
generalized  BFG  using  the  probabilistic  data  association  (PDA)  algorithm.  The  PDA  algorithm  starts  by 
assigning  prior  model  probabilities  (typically  equal)  for  all  the  possible  models  in  the  library  and 
recursively  updates  those  probabilities  based  on  the  voltage  and  current  measurements.  In  the  event  of 
an  unknown  battery  to  be  gauged,  the  PDA  algorithm  selects  the  most  similar  OCV  model  to  the  battery 
from  the  library.  We  also  demonstrate  a  strategy  to  select  the  minimum  sets  of  OCV  parameters  rep¬ 
resenting  a  large  number  of  Li-ion  batteries.  The  proposed  approaches  are  demonstrated  using  data  from 
portable  Li-ion  batteries. 


©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Rechargeable  batteries  have  increasingly  been  finding  applica¬ 
tions  in  many  emerging  frontiers;  consumer  electronics,  aerospace 
equipment,  electric  vehicles  and  smart  grid  are  some  examples. 
Lithium-ion  and  lithium-polymer  based  rechargeable  batteries  are 
widely  adopted  due  to  their  high  energy  density  and  long  life. 
Battery  fuel  gauging  (BFG)  plays  an  important  role  in  rechargeable 
battery  applications.  A  BFG  tracks  the  state  of  charge  (SOC)  and 
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state  of  health  (SOH)  of  the  battery  and  predicts  the  time  to  shut 
down  (TTS)  and  the  remaining  useful  life  (RUL). 

Li-ion/polymer  batteries  are  being  produced  by  numerous 
manufacturers  and  the  exact  chemistry  of  these  batteries  is  usually 
not  available  (mostly)  due  to  proprietary  reasons.  In  the  current 
business  practice,  battery  manufacturers  do  not  provide  BFG  ca¬ 
pabilities  —  it  is  up  to  the  device  manufacturer  to  develop  the  BFG 
which,  in  the  present  state  of  the  art,  requires  some  parameters  that 
are  battery  specific.  Further,  owing  to  the  fact  that  the  battery 
characteristics  drift  with  environmental  and  usage  factors,  a  BFG 
should  have  the  ability  to  perform  independently  of  the  manufac¬ 
turer  data-sheets.  A  robust  BFG  capability  allows  an  end  user  to 
rapidly  switch  batteries  without  requiring  the  user  to  adhere  to  a 
particular  battery  manufacturer.  This  is  particularly  useful  in 
portable  electronic  applications  where  battery  swaps  are  frequent. 
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Nomenclature 

a[/<]  battery  equivalent  ( voltage  drop )  model  (see  (2)) 

b  battery  equivalent  model  parameter  (2) 

BFG  battery  fuel  gauge 

ch  Coulomb  counting  coefficient  (1 ) 

FP[k]  derivative  of  OCV  a  (6) 

k j  OCV  parameters  of  thejth  model  (2) 

nd[k]  measurement  noise  of  the  voltage  drop  model  (2) 

OCV  open  circuit  voltage 

P r(*s[/<])  OCV  model  (2) 

Q.  process  noise  variance  (1 ) 

R  measurement  noise  variance  (of  the  voltage  drop 

model)  (2) 

d[k]  likelihood  of  jth  model  at  time  k  (14) 

SOC  state  of  charge 

ws[k]  process  noise  (1) 

Wj[k]  Kalman  filter  gain  corresponding  to  model  j  (18) 
xs[k]  state  of  charge  (SOC)  of  the  battery  (1 ) 

Zf[k\  measured  current  through  the  battery  (1 ) 

Zy[l<]  measured  terminal  voltage  (2) 

ft\k]  probability  of  thejth  model  at  time  /c  (13) 
v>[k]  Kalman  filter  innovation  corresponding  to  model  j 
(18) 


Most  of  the  existing  BFG  algorithms  are  based  on  an  OCV  model 
and  a  dynamic  equivalent  circuit  model  (ECM);  the  parameters  of 
the  OCV  curve  (shown  to  be  unchanged  over  temperature  changes 
and  aging  [1  )  are  estimated  offline  and  the  ECM  parameters 
(known  to  vary  with  temperature,  loading  [2]  and  aging)  are  esti¬ 
mated  online.  In  most  works,  (the  combined  model  [3,4  )  was 
adopted  for  modeling  the  OCV-SOC  relationship;  many  other 
possible  OCV  models  are  discussed  in  Ref.  [1  as  well.  In  Refs.  [3-8], 
different  dynamic  models  such  as  simple  (resistance  only)  model, 
zero/one  state  hysteresis  mode  and  enhanced  self  correcting 
models  are  discussed;  these  works  give  significant  attention  to 
modeling  the  hysteresis  effect  in  the  battery.  Combinations  of 
resistance/capacitance  models  where  the  hysteresis  effect  is  either 
ignored  or  modeled  as  an  error  can  be  found  in  Refs.  [9-18  . 

Based  on  the  application,  BFG  can  be  divided  into  two  cate¬ 
gories:  those  used  in  fixed  applications  and  those  used  in  portable 
applications.  In  fixed  applications,  the  BFG  has  to  continuously 
monitor  only  one  battery  (or  pack),  whereas  in  portable  applica¬ 
tions  the  BFG  might  be  expected  to  monitor  two  (or  more)  batteries 
—  one  at  a  time.  Table  1  compares  the  requirements  of  a  BFG  in 
fixed  vs.  mobile  applications. 

Our  focus  in  this  paper  is  on  portable  applications,  although  a 
generalized  BFG  will  be  useful  in  fixed  applications  as  well. 


Table  1 

Comparison  of  BFG  requirements  in  fixed  vs.  mobile  applications. 


Fixed  applications  Portable  applications 


Load  is  constantly 

connected  to  the  battery 

Battery  is  replaced  after  end  of  its  life; 
BFG  doesn't  have  to  keep  the 
parameters  of  the  old  battery. 


Examples:  electric  vehicles  (EV), 
power  tools,  aerospace  equipment. 


Possibility  for  battery  swap, 
e.g.,  a  smart  phone  user  keeping 
two  batteries  —  one  in  use  and 
one  as  back-up. 

Battery  swap  is  practiced,  i.e.,  the 
BFG  needs  to  keep  the  parameters 
of  two  (or  more)  batteries. 

When  a  back-up  battery  is  inserted, 
the  BFG  needs  to  adapt  to  the 
parameters  as  well. 

Examples:  smart  phones,  wearable 
devices,  tablets,  laptops,  etc. 


An  accurate  BFG  requires  the  offline  estimation  of  several  bat¬ 
tery  equivalent  model  parameters.  Recently,  we  developed  robust 
approaches  for  battery  fuel  gauging  that  require  the  characteriza¬ 
tion  of  only  the  open  circuit  voltage  (OCV)  of  the  battery  SOC 
[19,2,20  .  In  Refs.  [1  ,  the  many  aspects  of  OCV  parameter  estima¬ 
tion  of  Li-ion  battery  cells  are  elaborated.  Later  in  this  paper,  we 
show  that  the  OCV  parameters  of  several  types  of  batteries  can  be 
reduced  to  few  sets  of  OCV  parameters  though  machine  learning 
techniques. 

There  are  two  possible  approaches  to  develop  a  generalized  BFG 
that  is  independent  of  battery  chemistry.  The  first  approach  is  to 
simply  have  a  library  of  OCV  parameters  and  to  employ  online 
detection  approaches  in  order  to  select  the  most  suitable  OCV 
model  for  fuel  gauging.  In  other  words,  this  first  approach  seeks  to 
resolve  the  association  ambiguity  between  several  possible  OCV 
parameters  and  the  battery  being  monitored.  The  second  approach 
seeks  to  estimate  the  OCV  parameters  from  online  data  [21-23];  an 
iterative  process  is  employed  to  continuously  update  the  OCV  pa¬ 
rameters  and  battery  capacity.  One  of  the  challenges  faced  by  the 
second  approach  is  the  detection  of  new  batteries  (due  to  battery 
swaps)  and  restarting  the  OCV  parameter  estimation  routine  only 
when  required.  Further,  the  iterative  estimation  of  SOC,  OCV  &  ECM 
parameters  and  the  battery  capacity  adds  to  the  instability  of  the 
BFG  algorithm  and  costs  in  terms  of  robustness. 

The  present  paper  takes  the  first  approach  to  chemistry  adaptive 
battery  fuel  gauging.  Currently,  we  are  investigating  the  second 
approach  as  well.  To  the  knowledge  of  the  authors,  no  other  work 
has  attempted  to  develop  a  generalized  BFG  that  adapts  to  battery 
chemistries  online.  The  contributions  of  this  paper  are  listed  below: 

1.  Derivation  of  a  generalized  BFG  based  on  PDA.  We  derive  a  novel 
approach  to  chemistry-adaptive  battery  fuel  gauging  based  on 
the  principles  of  probabilistic  data  association  (24  .  As  the  name 
suggests,  the  PDA  assesses  the  probability  of  each  known  OCV 
model  against  the  observations  (instantaneous  voltage  and 
current)  and  employs  Gaussian  mixture  reduction  technique  to 
recursively  propagate  the  best  possible  posterior  probability 
density  of  the  SOC. 

2.  The  PDA  approach  also  can  serve  as  a  means  to  detect  a  new 
battery.  The  association  probability  can  also  serve  as  a  means  to 
new  battery  detection,  as  far  as  the  new  battery  significantly 
differs  from  the  current  battery  in  its  OCV  characteristics,  ca¬ 
pacity  and  resistance. 

3.  The  PDA  algorithm  enables  the  BFG  to  work  on  unknown  batteries. 
In  the  event  that  an  unknown  battery  for  which  the  OCV  pa¬ 
rameters  do  not  belong  to  any  of  the  ones  in  the  library,  the  PDA 
algorithm  will  find  the  most  similar  OCV  model  from  the  library 
and  use  it  to  gauge  the  unknown  battery;  assuming  that  the 
library  has  a  diverse  set  of  OCV  parameters,  the  PDA  based  BFG 
becomes  applicable  on  unknown  batteries. 

4.  A  clustering  approach  to  store  minimal  set  of  OCV  parameters  to  be 
used  by  the  generalized  BFG.  By  exploiting  machine  learning 
techniques,  we  demonstrate  a  novel  approach  to  select  the 
minimal  number  of  OCV  models  without  compromising  the  SOC 
tracking  performance. 

The  initial  interest  in  developing  universal  battery  solutions 
came  from  battery  charger  designers  in  order  to  address  the  ever 
increasing  number  of  different  types  of  battery  operated  equipment 
and  each  of  these  having  its  own  customized  charger;  this  increases 
the  amount  of  electronic  wastage  and  adds  to  the  cost  of  the  device. 
Hence,  there  has  been  significant  effort  in  developing  a  universal 


1  The  second,  online  OCV  estimation  approach  is  left  for  future  investigation. 
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battery  charger  [25-28  .  In  the  most  simplest  form,  the  first  chal¬ 
lenge  in  designing  a  universal  charger  is  to  find  the  voltage  to  which 
the  battery  has  to  be  charged.  The  approaches  in  Refs.  [25-28] 
suggest  using  a  look-up  table  of  incremental  voltage  in  response 
to  charging  by  a  certain  number  of  Coulombs. 

It  must  be  noted  that  the  universal  charging  approaches  dis¬ 
cussed  in  Refs.  [25-28]  attempt  to  charge  the  battery  in  a  very 
simple  form  -  finding  the  correct  voltage  to  apply.  The  optimal 
battery  charging  is  an  entirely  different  problem  (see  Refs.  [29-31  ] 
and  the  references  therein)  which  aims  to  apply  optimal  charging 
waveform  to  the  battery  such  that  the  charging  time  can  be  reduced 
without  affecting  the  longevity  of  the  battery.  The  optimal  charging 
algorithm  (OCA)  also  requires  accurate  knowledge  of  the  battery 
SOC  which  makes  BFG  capability  a  requirement  in  OCAs.  Our  focus 
on  this  paper  is  limited  to  developing  a  chemistry  adaptive  BFG. 

The  rest  of  this  paper  is  organized  as  follows:  In  section  2,  we 
develop  the  PDA  algorithm  for  the  robust  multimodel  BFG  with 
many  possible  OCV  model  parameters.  In  section  3,  we  show  how 
to  select  the  right  OCV  parameter  set  for  storage  from  among  many 
possible  candidates.  Computer  simulations  of  the  proposed  ap¬ 
proaches  are  presented  in  section  4  based  on  simulated  data  as  well 
as  real  battery  data  and  the  paper  is  concluded  in  section  5. 


zv  [k]  ~  VJ0(xs[k\k  -  1])  +  H>  [k]  (xs[k]  -  xs[k\k  -  1])  +  aT  [k\ b 
+  nd  [k]  =  [k]xs  [k]  +  e>  [k]  +  ar  [k]  b  +  nd  [k] 

where 

H>  [k]  = 

L  J  dx 


(7) 

x=xs[k\k-Y] 


e>  [k]  =  vl(xs[k\k  -  1])  -  Hj  [k]xs  [k\k  -  1 


Given  the  measured  voltage  zv[k]  and  current  z/[/<],  the  objective 
is  to  recursively  estimate  the  SOC  xs[k]. 

Using  Bayes  rule,  the  posterior  probability  density  of  xs[k]  given 
all  the  measurements  up  to  k,  zv[l:/<],  is  written  in  (9) 


Nj 

p(xs[k]\zv[\  :  k] )  =Y2p(Xs[k]\j,Zv[  1  :  k])p(j\Z„[\  :fc]) 

J=1 

_  Ej!iP(Z,M|Xs[k]J,Z„[l  :  k-  1])P0>,[1  :  k])p(Xs{k]\Zv[\  :  k-  1]) 


p(z„[k]|z„[l  :  k  —  1  ]  j) 


(9) 


2.  PDA  approach  for  robust  multimodel  BFG 


where  the  denominator  is  the  normalizing  factor  and 


The  SOC  tracking  problem  [19,2,20]  with  multiple  models  for 
the  OCV  parameter  set  is  summarized  below: 

xs[k]  =  xs[k  -  1]  +  chAZ([k  -  1]  +  ws[k]  (1) 


p(z„[k]|xs[k]  J,z„[l  :  k-  1]) 

=  J\f(zv  [k] ;  [k]  xs  [k]  +  ei  [k]  +  aT  [k]  b,  R 


(10) 


V]0(xs[k\)  +  aT  [k]b  +  nd  [k] 

p  (j j  z„  [1  k] )  =  0  [k] 

(11) 

pT (xs [fc] ) kj  +  aT [fc] b  +  nd[k]  j  =  1 , . . . , Nj 

(2)  p(xs[k]|z/l  :  fc- 1])  =Af(xs[k];xs[k|k-  l],ps[k|k-  1]) 

(12) 

where  xs[k]  is  the  SOC  of  the  battery,  Ch  is  the  Coulomb  counting 
coefficient  [2  ,  A  is  the  sampling  time,  zv[k]  is  the  measured  voltage, 
Zi[1<]  is  the  measured  current,  ws[k]  and  rid[k]  are  the  process  and 
measurement  noises  (see  Ref.  [20]  for  description),  which  are 
assumed  to  be  zero-mean  white  with  standard  deviations  ^/Q  and 
VR,  respectively,  b  corresponds  to  the  dynamic  equivalent  circuit 
model  parameter  vector  (which  has  model  dependent  dimensions, 
see  Ref.  19]),  a  is  a  known  vector  depending  on  the  dynamic 
equivalent  circuit  model  used,  and  the  OCV  model  is  given  by 


pr(xs[k]) 


1111 

(xs[k])2  (xs[k])3  (xs[k])4 * * 


The  model  probability^  /[fc]  is  updated  to  its  current  value  as. 


0  [k]  =p(j\zv[1:k])  = 


P(Zv[k]\j,Zv[1  :k-1])p(j\zv[1  :  k  —  1  ] ) 


p(Z„[k]|Z„[l  : fc —  1]) 

p(z„  [k]  ]j,Zv  [1  :  k  -  1  ]  )p(/|z„  [1  :  k  -  1  ] ) 
YfilP(Zv{k]\j,Zv{1  :k-1])p(j\Zv[1:k-1]) 

Cj[k]0[k- 1; 


Y^^[k]0[k- 1 


(13) 


where 


xs[k]  log(xs[k])  log(l  —  xs[k]) 


(3) 


d[k]  =p(z„[k][/,z„[  1  :  k 


and  the/11  set  of  OCV  model  parameters  are  given  by 


=  n(zv  [k] ;  2  [fc|k  -  1] ,  9  [fc] 

=  A f(o-y[k],9[k]) 


(14) 


and 


kje{ki>k2,  ...,kNj}  (5) 

indicates  Nj  OCV  model  parameters. 

The  measurement  model  (2)  can  be  written,  to  a  first  order 
approximation,  as 


2  Fast  charging  ages  the  battery  fast. 


is  the  likelihood  of  the  OCV  model  j  (specified  by  lq),  and 


2[k]  =  Zv  [fc]  -4[k|k-  1] 


S'  [fc]  =  (h>  [k] )  Ps  [k|k  -  1]  +  R 


(15) 

(16) 


3  In  the  PDAF  for  tracking  in  clutter  [24]  this  is  the  (association)  probability  of 

measurement  j  originating  from  the  target  of  interest;  here,  this  is  the  probability 

that  the  (only)  measurement  originated  form  model  j.  Consequently,  we  use  the 

designation  of  model  probability. 
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are  the  innovation  and  the  variance  of  the  innovation, 
respectively. 

Now,  (9)  is  explicitly  written  as 


Ni 


p (xs  [k]  \zv[\  :  k] )  [k] Af  (zv [k]  ■  PP  [k] xs  [fc]  +  0  [k] , R 

j=  i 


•A/'^xs[k];xs[k|k -  l],Ps[k|k -  1]) 


Ni 


=  c  Y)  0  [/<]  AT  ( xs  [/<] ;  x1  [k  |  /<] ,  ^  [k  |  k] 
j=  i 

where  c  is  the  normalizing  constant  and 

xj.  [k | k\  =  xs  [k\k  -  1]  +  W-i  [k]  ft  [k\ 

f^[k|k]  =  Ps[k\k-  1]  -  (wj [k])2 S' [k] 

=  (l  -w,[k]H,[k])  Ps[k\k-1]  +  wj\k]2R 


(18) 


are  the  “mode  matched”  SOC  estimates  and  the  corresponding 
estimation  error  variances,  and 


Wj[k }  = 


Ps[k\k-  l]HJ[k] 

W\ 


is  the  Kalman  gain  corresponding  to  thejth  model. 

The  posterior  distribution  p(xs[/c]|zv[l  :/<])  becomes  a  Gaussian 
mixture  as  shown  in  (17).  In  recursive  estimation,  we  seek  to 
approximate  this  mixture  as  a  single  Gaussian  component.  There 
are  several  approaches  for  Gaussian  mixture  reduction  in  the  liter¬ 
ature,  such  as  [32]  and  the  references  therein.  It  is  shown  in  Refs. 
[[33],Theorem2]  that,  the  Gaussian  density  Af{xs[k\;xs[k\k\,Ps[k\k]) 
minimizes  the  Kullback-Leibler  distance  between  (17)  and  its 
Gaussian  approximation,  where 


NJ  .  NJ 

xs  [k | k\  =  ^  ft  [k] ^  [k | k\  =  xs  [k\k  -  1]  +  ^  ft  [k]  WJ  [k\  ft  [k] 

j= i  j=  i 


(21) 


Ps[k\k] 


Nj 

Y)  0  [k]  fpi  [k|k]  +  (p4  [k|k]  -  xs  [k|k] 
j=i 

nj  .  r  ■  •  /  • 

52 0 [k]  P>s [k\k]  +  w’ [k] (v> [k]  -0k] 

j= i  L 


(22) 


where 


Then,  the  following  mode-matched  filtering  and  likelihood 
computations  are  performed  for.  j  =  1,2,...N/: 


Meas.  prediction  z1 


k 


k- 1 


-  pT  ( xs  [k  |  k  - 1  ] )  kf  +  aT  [k\  b  (26) 


Innovation  ft  [k]  =  z[k\  -z1 


Linearized  H  model  HJ  [/<] 


k 


k-  1 


(27) 


dVj(x) 

dx 


x=xs  [k\k-\] 


(^7)  Innovation  variance  ft  [k\  =  R  +  [k\ )  Ps  [k\k  -  1" 


Filter  gain  WJ  [k]  = 


Ps[k\k-\]Hj[k } 
R+  (HJ  [k])2 Ps[k\k  -  1 


(28) 


(29) 


(30) 


/  -j  g  \  Estimation  variance  0  \k\k]  -  ( 1  -  W’  \  k  HJ  [  k] )  Ps[k|/c-l 


+  (wfk] )  R 


(31) 


(20)  Model  likelihood  0[k]  =  Af (z [k];0[k\k-  1], 0 [/<])  (32) 


Then,  the  model  probabilities  are  updated  as 


0  [k]  1] 


(33) 


Finally,  the  SOC  estimate  and  the  estimation  error  variance  are 
given  by 


xs  [k\k]  =  xs  [k\k  -  1]  +  ^  Wj  [/<]  ft  [k\  ft  [k] 

j=  i 

Nj  r  2 

Ps[k|k]  =X)/?[k]  pj[k|k]  +  [l<]  [l<] -v[k]) 

j= i  L 


(34) 


(35) 


Remark :  It  must  be  noted  that  the  above  procedure  is  a  novel  use 
of  the  probabilistic  data  association  (PDA)  filter  from  Ref.  24]  for 
the  purpose  of  model  identification. 

A  block  diagram  of  the  proposed  approach  is  shown  in  Fig.  1. 

3.  Machine  learning  approach  to  efficient  OCV  parameter 
storage 


v  [k]  =  Y  0  M  wj  M  j  M  (23) 

i=i 

Now,  starting  from  the  SOC  estimate  xs[k  —  1 1  /<  —  1  ]  and  the 
corresponding  estimation  error  variance  Ps[/<-l|/c-l],  the  updated 
SOC  estimate  xs[k\k]  and  the  corresponding  estimation  error  vari¬ 
ance  Ps[/c|/<]  can  be  obtained  as  follows  [34]: 

First,  the  state  prediction  and  the  state  prediction  variance,  both 
of  which  do  not  depend  on  the  OCV  models,  are  computed  as. 

Prediction  xs[k\k  -  1]  =  xs[k  -  1  \k  -  1]  +  chAZi[k]  (24) 

Prediction  cov.  Ps[k\k  -  1]  =  P[k  -  1  \k  -  1]  +  Q.  (25) 


BFGs  in  portable  applications  are  preferred  to  have  the  ability 
gauge  any  battery  without  strictly  requiring  offline  OCV  charac¬ 
terization.  This  is  partly  due  to  the  time  and  cost  involved  in  the 
OCV  characterization  process;  the  possibility  of  frequent  battery 
changes  (including  swaps)  also  make  such  BFGs  very  desirable. 

It  is  generally  understood  that  the  shape  of  the  OCV  curve  de¬ 
pends  on  the  chemical  and  physical  structure  of  the  battery. 
However,  there  are  no  analytical  expressions  for  the  OCV  charac¬ 
terization  in  terms  of  chemical  and  physical  structure.  In  practice, 
the  OCV  curves  are  obtained  through  offline  experimentation.  In 
Ref.  [1]  we  review  many  existing  approaches  for  OCV  parameter 
estimation  trough  experiment.  It  is  found  that  the  OCV  curves 
remain  unchanged  across  temperature  and  that  the  OCV  curves  of 
many  batteries  by  different  manufacturers  are  similar  (provided 
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that  the  normalized  OCV  modeling  approach  described  in  Ref.  [1]  is 
followed).  Hence,  it  must  be  possible  for  a  BFG  to  store  a  few  sets  of 
OCV  parameters  and  be  able  to  gauge  batteries  made  by  many 
different  manufacturers. 

In  this  section,  we  discuss  ways  to  select  the  OCV  parameters 
{k/ljl-!  that  are  used  in  the  multimodel  BFG.  Table  2  shows  the 
details  of  the  battery  cells  that  were  used  in  the  experiment;  the 
fourth  column  titled  “Cell  Number”  indicates  the  different  cells  of 
the  same  battery  type  used  in  the  experiment. 

The  OCV  characterization  experiment  1]  involves  the  following 
four  sequential  steps: 

1.  Fully  charging  the  battery  at  room  temperature 

2.  Bringing  the  battery  to  test  temperature 

3.  Fully  dischargnig  the  battery  at  a  slow  rate  (usually  a  C/30  rate) 
at  the  test  temperature,  and 

4.  Charging  the  battery  back,  until  full,  at  a  slow  rate  (usually  a  C/ 
30  rate)  at  the  same  test  temperature 

The  approach  described  in  Ref.  [1  can  be  employed  to  estimate 
the  OCV  parameters.  The  OCV  characterization  in  the  above  manner 
was  performed  at  15  different  temperatures  from  Ref.  -25  °C  to 


45  °C  in  steps  of  5  °C.  Hence,  a  total  of  34  x  15  =  510  OCV  charac¬ 
terization  experiments  were  performed.  In  Refs.  [1  ,  we  provide  a 
complete  analysis  of  these  OCV  experiments. 

One  of  the  important  findings  from  Ref.  1  is  that  the  OCV-SOC 
characterization  curve  can  be  made  to  be  the  same  across  tem¬ 
peratures  by  using  the  normalized  OCV  modeling.  In  the  normalized 
OCV  modeling,  the  battery  capacity  is  computed  at  each  tempera¬ 
ture  and  used  for  the  computation  of  SOC  resulting  in  a  single 
OCV-SOC  curve  across  all  temperatures.  This  enables  significant 
savings  in  OCV  parameter  storage,  i.e.,  it  is  sufficient  to  store  a 
single  set  of  OCV  parameters  for  each  battery  type  and  estimate  the 
battery  capacity  online  [2]. 

In  Table  2,  we  show  the  estimated  OCV  parameters  for  each  cell. 
These  parameters  can  be  used  to  generate  the  OCV  for  a  particular 
SOC  as 

Vo(s)  =  P(s)rk  (36) 

Fig.  2  shows  the  OCV-SOC  curves  of  cell  listed  in  able  2.  It  can 
be  noticed  that,  regardless  of  the  battery  type,  many  of  the  OCV 
curves  are  similar.  Indeed,  for  many  battery  types,  the  differences 
between  their  OCV  curves  are  smaller  than  the  error  in  computing 


Measurement 

zv[k] 


Fig.  1.  Block  diagram  of  the  proposed  multimodel  BFG. 
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Table  2 

List  of  Batteries  and  their  OCV  Parameters. 


Sample  no. 

Battery  manufacturer 

Serial  number 

Cell  number 

K0 

K1 

K2 

I<3 

K4 

K5 

K6 

K7 

K8 

1 

Nokia 

BP  4L 

1 

-4.4813 

30.3625 

-5.4895 

0.6087 

-0.0285 

-16.9776 

38.2121 

-0.1162 

0.5813 

2 

Nokia 

BP  4L 

2 

-4.4882 

34.7351 

-6.1389 

0.668 

-0.0308 

-20.8428 

44.9474 

-0.1444 

0.5581 

3 

Nokia 

BP  4L 

3 

-5.7439 

31.3426 

-5.7249 

0.6379 

-0.0299 

-16.2129 

38.615 

-0.0134 

0.5889 

4 

Nokia 

BP  4L 

4 

-4.1289 

45.3504 

-7.7859 

0.8302 

-0.0376 

-30.8211 

61.2808 

-0.3192 

0.3708 

5 

Samsung 

EB555157VA 

1 

-1.5085 

41.2402 

-7.07 

0.7681 

-0.0359 

-30.3721 

57.0893 

-0.4572 

0.2944 

6 

Samsung 

EB555157VA 

2 

-2.2216 

38.8595 

-6.7447 

0.7396 

-0.0348 

-27.3031 

52.8784 

-0.3599 

0.3395 

7 

Samsung 

AB463651 

1 

-2.7886 

36.6116 

-6.4458 

0.7156 

-0.0341 

-24.5508 

48.9357 

-0.2698 

0.4056 

8 

Samsung 

AB463651 

2 

-1.9049 

34.9858 

-6.1037 

0.6732 

-0.0319 

-24.2116 

47.3439 

-0.316 

0.492 

9 

Samsung 

AB463651 

3 

-2.3237 

47.4112 

-8.1229 

0.8811 

-0.0411 

-34.8322 

65.6118 

-0.4617 

0.4508 

10 

Samsung 

AB463651 

4 

-1.3313 

28.2101 

-4.8784 

0.5292 

-0.0247 

-18.9233 

37.9843 

-0.27 

0.4251 

11 

Samsung 

EB504465 

1 

-2.3739 

32.6716 

-5.784 

0.6456 

-0.0309 

-21.5345 

43.4031 

-0.2538 

0.3799 

12 

Samsung 

EB504465 

2 

-1.5684 

46.4157 

-7.9021 

0.8536 

-0.0397 

-34.918 

64.8046 

-0.5132 

0.3358 

13 

Samsung 

EB504465 

3 

-1.962 

33.7161 

-5.917 

0.656 

-0.0313 

-23.0268 

45.3406 

-0.3099 

0.3483 

14 

Samsung 

EB504465 

4 

-1.6756 

42.3832 

-7.274 

0.7917 

-0.0371 

-31.1357 

58.6333 

-0.4389 

0.356 

15 

Samsung 

EB575152 

1 

-2.0542 

28.6842 

-4.9746 

0.5372 

-0.0249 

-18.3279 

38.1599 

-0.1554 

0.6502 

16 

Samsung 

EB575152 

2 

-1.9791 

36.6886 

-6.2418 

0.6632 

-0.0302 

-25.6191 

50.2428 

-0.2748 

0.6608 

17 

Samsung 

EB575152 

3 

-2.5887 

36.6903 

-6.2908 

0.6715 

-0.0306 

-24.8732 

49.677 

-0.2402 

0.6835 

18 

Samsung 

EB575152 

4 

-2.9254 

25.4594 

-4.6864 

0.539 

-0.0265 

-14.3882 

32.0363 

-0.1311 

0.3392 

19 

Blackerry 

FS  1 

1 

-2.868 

19.5393 

-3.7414 

0.4441 

-0.0224 

-9.1018 

23.1865 

-0.0319 

0.4245 

20 

Blackerry 

FS  1 

2 

-1.6544 

33.2971 

-5.8206 

0.6441 

-0.0307 

-23.0123 

45.0151 

-0.3222 

0.4157 

21 

Blackerry 

M  SI 

1 

-2.2305 

24.9583 

-4.4232 

0.4878 

-0.0231 

-14.8639 

32.3114 

-0.1573 

0.438 

22 

Blackerry 

M  SI 

2 

-2.1679 

25.4322 

-4.4909 

0.4933 

-0.0232 

-15.4082 

33.0761 

-0.1848 

0.4029 

23 

Blackerry 

M  SI 

3 

-2.1511 

28.6238 

-4.9981 

0.5446 

-0.0255 

-18.2572 

37.8914 

-0.2113 

0.4905 

24 

Blackerry 

M  SI 

4 

-0.9741 

30.3903 

-5.2017 

0.5596 

-0.0259 

-21.3917 

41.5533 

-0.3517 

0.3988 

25 

LG 

LGIP 

1 

-3.3584 

36.5011 

-6.4703 

0.7216 

-0.0344 

-23.6683 

48.2804 

-0.214 

0.4695 

26 

LG 

LGIP 

3 

-2.4388 

38.5035 

-6.718 

0.7406 

-0.035 

-26.6555 

52.1114 

-0.3184 

0.4415 

27 

Huawei 

HB4Q1 

1 

-1.6226 

39.0734 

-6.722 

0.732 

-0.0343 

-28.27 

53.7686 

-0.4079 

0.3679 

28 

Huawei 

HB4Q1 

2 

-2.4954 

35.1997 

-6.1698 

0.6815 

-0.0323 

-23.6582 

47.1783 

-0.2781 

0.395 

29 

Samsung 

EBL1A2GBA 

1 

-1.9747 

50.2979 

-8.5394 

0.9183 

-0.0424 

-37.858 

70.3186 

-0.522 

0.3504 

30 

Samsung 

EBL1A2GBA 

2 

-1.8145 

50.823 

-8.6133 

0.9251 

-0.0427 

-38.5429 

71.2411 

-0.5464 

0.3281 

31 

Samsung 

EBL1A2GBA 

3 

-1.611 

49.5084 

-8.3989 

0.9035 

-0.0418 

-37.6603 

69.4047 

-0.5627 

0.2877 

32 

Samsung 

EBL1A2GBA 

4 

-2.276 

49.3495 

-8.4102 

0.907 

-0.042 

-36.6127 

68.6301 

-0.4817 

0.3443 

33 

Samsung 

L1G6LLAGS 

1 

-4.4256 

83.6175 

-14.0185 

1.4878 

-0.0678 

-64.5844 

118.3224 

-0.8474 

0.3087 

34 

Samsung 

L1G6LLAGS 

2 

-4.1835 

83.7601 

-14.0218 

1.4864 

-0.0676 

-65.0311 

118.7525 

-0.871 

0.2839 

these  curves  [1  .  Hence,  the  storage  of  OCV  parameters  can  be 
further  reduced  by  selecting  curves  that  represent  many  similar 
batteries. 

A  simple  approach  to  reduce  the  storage  of  OCV  parameters 
involves  the  following  steps: 


groups  found  based  on  the  analysis  of  parameters  listed  in 
Table  2.  These  three  groups  are  labeled  as  follows:  OCV1  - 
green;  OCV2  -  red;  and  OCV3  -  blue. 

The  above  approach  is  summarized  in  Fig.  5. 


1.  Collect  the  OCV  parameters  from  as  many  different  battery  cells  as 
possible :  It  is  important  to  collect  the  OCV  parameters  so  that 
they  span  the  envelope  of  possible  OCV-SOC  curves.  This  helps 
in  reducing  errors  in  adaptive  BFG.  Multiple  cells  of  the  same 
battery  in  Table  2  serve  this  purpose. 

2.  Cluster  the  OCV  parameters  based  on  similarity :  The  similarity  of 
OCV  parameters  can  be  simply  computed  be  comparing  the  OCV 
parameters  of  each  cell.  We  compute  the  Euclidean  distance  as  a 
measure  of  dissimilarity,  i.e.,  the  dissimilarity  between  k/  and  k j 
is  given  by 


dij=  £  (Pr(s)ki-PT(s)kj)2  (37) 

V  set°T] 

Fig.  3  shows  the  dissimilarities  of  parameters  in  fable  2  as  a 
dendrogram  which  is  an  output  of  hierarchical  clustering  [35]. 

3  Select  one  OCV  parameter  for  each  cluster :  Once  the  OCV  pa¬ 
rameters  are  grouped,  we  take  the  average  in  each  group  and 
store  them  as  the  OCV  parameter  representing  that  group.  This 
particular  parameter  is  specified  as  k j  corresponding  to  the  jth 
group  as  used  in  the  multimodel  BFG.  Fig.  4  shows  the  three 


4.  Simulation  results 

In  this  section,  we  demonstrate  the  performance  of  the  pro¬ 
posed  algorithm  on  actual  battery  data.  The  multimodel  BFG  is 


4 


Maximum  likelihood  combination  of  OCV  parameters  is  an  alternative  method. 


Fig.  2.  OCV  curves  of  each  battery  cell  listed  in  Table  2. 
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Fig.  3.  Hierarchical  clustering  of  OCV  curves.  Dissimilarity  among  the  OCV  parameters 
are  displayed  as  a  dendrogram.  It  is  easy  to  see  three  groups  of  OCV  parameters. 

equipped  with  three  OCV  parameters  corresponding  to  the  three 
curves  shown  in  Fig.  4.  For  simplicity,  we  label  these  curves  as 
OCV1,  OCV2  and  OCV3.  Table  3  summarizes  the  details  of  these 
curves. 

The  multimodel  BFG  algorithm  proposed  in  this  paper  is 
implemented  through  state-control  in  some  applications,  such  as  in 
smart  phones,  the  BFG  needs  to  face  such  events  as  battery  swaps, 
new  battery  replacement  and  unconventional  usage,  such  as  leav¬ 
ing  the  battery  unused  until  it  becomes  “dead”;  using  non-standard 
chargers  that  may  apply  unusual  voltages;  etc.  The  state-control 
[36]  checks  the  voltage  and  current  data  for  unfavorable  condi¬ 
tions  (e.g.,  constant  current  is  not  favorable  for  parameter  estima¬ 
tion)  and  activates,  pauses,  resets  and  reactivates  different  modules 
of  the  BFG,  such  as  parameter  estimation  module,  capacity  esti¬ 
mation  module  and  SOC  tracking  module. 

Now,  we  demonstrate  the  performance  of  the  multimodel  BFG 
through  three  different  examples  in  the  next  three  subsections. 

4.2.  Demonstration  of  PDA  in  resolving  battery  type 

First,  we  demonstrate  the  ability  of  PDA  in  BFG  to  adapt  to 
different  battery  chemistries  without  the  exact  knowledge  of  the 
OCV  parameters.  Fig.  6  shows  the  data  used  for  this  experiment. 
This  data  was  collected  from  a  Samsung  battery  with  serial  number 
EB575152  as  follows:  Starting  from  a  fully  charged  battery,  dynamic 


Fig.  4.  Grouped  OCV  curves:  The  top  plot  shows  the  OCV  curves  corresponding  to  each 
group  in  Fig.  3  in  the  same  color.  The  plot  a  the  bottom  shows  the  representative  OCV 
plot  selected  for  each  group.  We  label  the  three  plots  as  follows:  OCV1  -  green;  OCV2 
-  red;  and  OCV3  -  blue.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


!VoV: 


Fig.  5.  Proposed  approach  for  efficient  parameter  storage. 
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Table  3 

Identification  of  three  OCV  curves. 


Label 

Max.  OCV 

Color/position  in  Fig.  4 

OCV1 

4.17  V 

Blue/lower 

OCV2 

4.17  V 

Red/middle 

OCV3 

4.35  V 

Green/top 

discharge  profiles  as  shown  in  Fig.  6  (bottom)  were  applied  for  6  h; 
After  that,  the  batter  was  charged  using  a  constant  current  of 
750  mA  for  30  min;  Finally,  the  dynamic  profile  was  applied  until 
the  voltage  across  the  battery  terminals  reached  the  shut  down 
voltage  of  3  V. 

The  voltage  and  current  data  shown  in  Fig.  6  were  fed  to  the 
multimodel  BFG  algorithm  in  order  to  track  the  battery  SOC.  The 
(multimodel)  BFG  algorithm  employs  three  sets  of  OCV  parameters 
listed  in  Table  3.  The  SOC  tracking  performance  of  the  multimodel 
BFG  is  shown  in  the  bottom  of  Fig.  7.  Initially,  the  PDA  algorithm  in 
the  multimodel  BFG  assigns  equal  probabilities  (1/3)  for  each  OCV 
model.  The  correct  OCV  model  corresponding  to  the  input  data  is 
identified  within  (approximately)  1  h  of  the  start  of  the  BFG  algo¬ 
rithm  —  the  mode  probability  of  OCV2  rises  to  1  in  the  top  plot  of 
Fig.  7. 

It  must  be  noted  that  the  BFG  algorithm  needs  an  initial  “guess” 
for  the  SOC  at  the  very  start.  It  can  be  observed  that  the  initial 
voltage  in  the  data  is  approximately  4.17  V  (see  the  top  plot  of 
Fig.  6.)  The  OCV  curves  OCV1  and  OCV2  correspond  to  100%  SOC  at 
4.17  V  and  OCV1  correspond  to  only  85%  SOC  at  4.17  V  (see  the 
bottom  plot  of  Fig.  4.)  The  initial  SOC  guess  then  becomes. 
85  +  100  +  100/3  =  95%. 

The  state-control  of  the  BFG  assumes  every  start  of  the  simu¬ 
lation  as  the  very  first  time  a  battery  is  seen  by  the  BFG.  Hence,  it 
takes  some  time  to  check  for  voltage/current  patterns  before  the 
parameter  estimation  and  PDA  are  activated.  Fig.  8  shows  a  close  up 
look  of  the  instance  when  the  PDA  algorithm  starts  its  computa¬ 
tions.  The  mode  probability  can  be  seen  updated  to  the  correct 
value  in  two  steps. 
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4.2.  Demonstration  of  PDA  in  portable  applications 

An  important  application  of  multimodel  BFG  arises  in  portable 
battery  applications  (see  Table  1).  Consider  smartphone  scenario 


Fig.  7.  Demonstration  of  multimodel  BFG  in  resolving  battery  type.  Initially,  the  BFG 
assigns  equal  probability  for  the  three  possible  OCV  models.  Due  to  the  similarity  of  the 
OCV1  and  OCV2  curves  at  high  SOC  region  (see  bottom  of  Fig.  4). 
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Fig.  6.  Voltage  and  current  from  the  battery.  The  data  is  collected  from  cell  4  of 
Samsung  EB575152  (see  Table  2). 


Time  (Hr) 


Fig.  8.  Mode  probability  update.  Markers  indicate  the  sampling  instances.  Once  the 
PDA  algorithm  starts,  the  accurate  mode  probability  is  reached  just  with  two  samples 
(0.2  s). 
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Fig.  9.  Battery  exchange  data  in  portable  applications. 


where  a  user  carries  a  fully  charged  second  (back-up)  battery  and 
when  the  current  battery  empties  it  is  replaced  by  the  back  up 
battery.  Fig.  9  shows  the  data  corresponding  to  the  above  battery 
replacement  scenario:  The  top  plot  shows  that  the  voltage  around 
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Fig.  11.  Estimated  parameters  in  the  replacement  scenario.  A  brief  dip  in  estimated 
parameters  can  be  observed  when  the  new  battery  is  inserted.  The  state-control  of  the 
BFG  routinely  checks  the  estimated  parameters  for  instabilities. 


the  14  h  mark  jumping  back  to  full  battery  voltage  with  no  charging 
current  present.  This  indicates  the  scenario  that,  once  the  current 
battery  becomes  empty,  it  is  replaced  by  a  different  battery  which 
might  not  have  the  same  type  of  OCV  curve. 

Fig.  10  shows  the  SOC  tracking  performance  of  the  multimodel 
BFG  and  the  corresponding  mode  probability.  Similar  to  the 
example  in  Section  4.1,  the  SOC  initialization  is  obtained  by 
combining  the  SOC-lookup  from  all  OCV  curves  in  the  multimodel 
BFG.  The  BFG  can  be  similarly  reinitialized  when  a  new  battery  is 
inserted;  however,  such  reinitialization  is  not  performed  in  Fig.  10. 
It  can  be  seen  that  the  BFG  is  able  to  adjust  to  the  correct  OCV  mode 
within  1  h  of  processing. 

The  BFG  is  implemented  with  the  assumption  of  a  dynamic 
equivalent  model  consisting  of  one  series  resistor  and  one  RC 
component;  this  would  result  in  a  4  x  1  model  parameter  vector  b 
(see  Ref.  19]  for  details).  Fig.  11  shows  the  dynamic  battery 
equivalent  model  parameter  estimates  for  the  example  corre¬ 
sponding  to  Fig.  10.  Each  element  of  b  is  labeled  as  b(l),  b(2),  b(3) 
and  b(4).  It  can  be  seen  that  the  new  battery  insertion  affects  the 
parameter  estimates  as  well. 


Fig.  10.  Demonstration  PDA  in  the  replacement  scenario. 


Fig.  12.  OCV  curves  used  in  simulated  PDA  experimentation. 
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Table  4 

Parameters  of  baseline  and  mixed  OCV  curves. 


Parameters 

ICO 

K1 

K2 

K3 

K4 

K5 

I<6 

K7 

OCV1 

-3.0927 

43.3102 

-7.4126 

0.7908 

-0.0359 

-30.3610 

59.0316 

-0.3932 

OCV2 

-1.2528 

49.7952 

-8.4537 

0.9127 

-0.0425 

-38.3972 

70.0370 

-0.6151 

OCV3 

-2.5215 

67.8629 

-11.3633 

1.2086 

-0.0554 

-52.9221 

96.3480 

-0.8136 

Mixed-OCV2 

-2.2678 

64.2493 

-10.7814 

1.1494 

-0.0528 

-50.0171 

91.0858 

-0.7739 

Mixed-OCV3 

-1.5065 

53.4087 

-9.0356 

0.9719 

-0.0451 

-41.3022 

75.2992 

-0.6548 

4.3.  Demonstration  of  PDA  when  the  battery  chemistry  is  unknown 

This  simulated  example  is  designed  to  demonstrate  the  ability  of 
the  PDA  algorithm  to  find  the  closest  OCV  parameter  set  from  the 
library  in  the  absence  of  exact  OCV  parameter  sets  corresponding  to 
the  battery  being  gauged.  The  mode  probability  in  the  PDA  algorithm 
will  converge  to  the  closest  OCV  parameter  set  that  gives  the  highest 
likelihood  for  the  data;  this  phenomenon  was  proved  in  Ref.  37]. 

In  Fig.  12,  we  introduce  two  more  OCV  curves  in  addition  to  the 
three  baseline  OCV  curves  (see  the  bottom  plot  in  Fig.  4)  that  we 
selected  through  the  machine  learning  approach  described  in 
Section  3;  parameters  for  these  two  additional  OCV  curves  are 
generated  as  follows. 

kmixed-OCV2  =  0.8kOCV2  +  0-2koCV3  (38) 

kmixed-OCV3  =  0.2kOCV2  +  0.8kOcV3  (39) 

where  each  k  is  an  OCV  parameter  set  corresponding  to  a  category 
indicated  in  the  subscript;  the  OCV  curve  corresponding  to  mixed- 
OCV2  is  more  similar  to  OCV2  than  it  is  to  OCV3;  and  the  OCV  curve 
corresponding  to  mixed-OCV3  is  more  similar  to  OCV3  than  it  is  to 
OCV2.  Table  4  lists  the  OCV  parameter  sets  corresponding  to  curves 
denoted  as  OCV1,  OCV2,  OCV3,  mixed-OCV2  and  mixed-OCV3  in 
Fig.  12. 

Now,  we  will  test  the  PDA-BFG  by  feeding  it  with  data  generated 
based  on  different  OCV  curves  shown  in  Fig.  12.  Detailed  de¬ 
scriptions  of  how  to  simulate  battery  data  are  discussed  in  Refs. 
[2,19,20  .  In  short,  first  the  true  SOC  is  computed  based  on  an 
assumed  battery  capacity  and  known  initial  SOC  (we  use  s[0]  =  1) 
for  a  dynamic  current  profile  similar  to  the  one  shown  in  Fig.  6; 


then,  the  observation  model  (2)  is  used  to  generate  voltage  mea¬ 
surements  with  an  assumed  voltage  measurement  error  s.d.  of 
Gy  —  1  mV  and  current  measurement  error  s.d.  of  ov  =  1  mA. 

First,  we  demonstrate  the  expected  SOC  tracking  performance 
when  the  exact  OCV  parameters  of  the  battery  being  gauged  is 
known  to  the  BFG.  Fig.  13  shows  the  summary  of  results  for  this 
scenario;  the  battery  data  was  simulated  using  the  parameters  of 
OCV3  curve  in  Fig.  12;  the  PDA-BFG  algorithm  has  the  OCV  pa¬ 
rameters  corresponding  to  curves  labeled  OCV1,  OCV2  and  OCV3  in 
Fig.  12.  The  top  plot  in  Fig.  13  shows  the  mode  probability  quickly 
converging  towards  the  correct  OCV  model;  the  second  and  third 
plots  show  the  SOC  tracking  error  quickly  converging  towards  zero. 

In  Fig  14,  the  SOC  tracking  performance  of  the  PDA-BFG  in  the 
absence  of  the  exact  OCV  parameter  is  demonstrated.  The  battery 
data  was  simulated  using  the  parameters  of  mixed-OCV2  curve  in 
Fig.  12;  as  before,  the  PDA-BFG  algorithm  has  only  the  OCV  pa¬ 
rameters  corresponding  to  curves  labeled  OCV1,  OCV2  and  OCV3  in 
Fig.  12.  The  parameters  of  mixed-OCV2  is  more  similar  to  OCV2 
than  any  other  OCV  parameters  in  the  library  (see  (38));  as  such,  the 
mode  probability  (top  plot)  of  PDA-BFG  quickly  converges  towards 
OCV2.  The  SOC  tracking  error  is  not  zero  (but  small  enough  -  less 
than  1%)  due  to  the  fact  that  the  OCV2  curve  is  close  but  not  the 
same  as  that  of  mixed-OCV2. 

Fig  15  shows  a  similar  mismatched  OCV  simulation  and  the 
corresponding  results  from  the  PDA-BFG  algorithm;  in  this  figure, 
the  data  is  generated  based  on  mixed-OCV3. 

5.  Conclusions 

We  presented  a  novel  approach  for  battery  fuel  gauging  that  is 
able  to  adapt  to  many  different  battery  chemistries.  We  showed 
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Fig.  13.  Demonstration  of  PDA-BFG  on  simulated  data.  The  data  is  generated  based  on  Fig.  14.  Demonstration  of  PDA-BFG  on  simulated  data.  The  data  is  generated  based  on  the 

the  parameters  of  OCV3  which  is  available  in  the  parameter  library  of  PDA-BFG.  parameters  of  mixed-OCV2  which  is  not  available  in  the  parameter  library  of  PDA-BFG. 
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Fig.  15.  Demonstration  of  PDA-BFG  on  simulated  data.  The  data  is  generated  based  on 
the  parameters  of  mixed-OCV3  which  is  not  available  in  the  parameter  library  of  PDA- 
BFG. 


in  a  recent  work  [1  ,  via  normalized  OCV  modeling,  that  the  OCV 
characterization  curve  does  not  change  with  the  temperature  of 
the  cell  -  which  enables  battery  fuel  gauging  (BFG)  with  a  single 
OCV  curve  regardless  of  the  temperature.  In  this  paper,  we  show 
that  the  (normalized)  OCV  curves  of  batteries  from  many 
different  manufacturers  are  similar.  Based  on  this  observation, 
and  by  using  machine  learning  algorithms,  we  develop  an 
approach  to  analyze  a  large  number  of  OCV  parameters  (repre¬ 
senting  many  different  batteries)  and  to  select  a  few  sets  of  pa¬ 
rameters  that  are  representative  of  the  entire  set  of  batteries.  The 
proposed  multimodel  BFG  approach  then  uses  those  selected 
OCV  parameters  and  the  probabilistic  data  association  (PDA)  al¬ 
gorithm  in  order  to  associate  the  battery  to  the  correct  OCV 
model,  i.e.,  identify  it.  The  resulting  multimodel  BFG  can  be  used 
to  gauge  a  wide  range  of  batteries  without  requiring  battery- 
specific  OCV  parameters. 
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Appendix  A.  Supplementary  data 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
dx.doi.org/10.1016/j.jpowsour.2014.09.006. 
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